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ABSTRACT 

Aims. Study the high-mass X-ray binary X-7 in M33 using broad-band optical data. 

Methods. We used recently published CFHT r' and i' data for variable stars in M33 to extract the light 
curve of the optical counterpart of X-7. We combined these data with DIRECT B and V measurements in 
order to search for an independent optical modulation with the X-ray periodicity. The periodic modulation 
is modelled with the ellipsoidal effect. We used UBVRr'i' magnitudes of the system to constrain the tem- 
perature and radius of the optical component. 

Results. The optical data revealed a periodicity of 3.4530 ± 0.0014 days , which is consistent with the 
known X-ray period. Double modulation, which we attributed to ellipsoidal modulation, is clearly seen in 
four different optical bands. The absolute magnitude in six optical bands is most consistent with a stellar 
counterpart with 33000 < T eff < 47000 K and 15 < R < 20%. We modelled the optical periodic modula- 
tion and derived the masses of the two components as a function of the orbital inclination and the radius of 
the stellar component. The resulting mass range for the compact object is 1.3 < M < 23 Mq. 
Conclusions. The system is probably a black hole HMXB, similar to Cyg X-l, LMC X-l and LMC X-3. 
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1. Introduction 

M33 X-7 (hereaft er X-7) is a well known bright variable X-ray source in th e nearb y M33 galaxy 



dLong et al 



1981I) . Its variability was identified as periodic by 



Peres et al 



dl989l) . who derived 



a period of 1.7857 d ays and were the firs t to su ggest its XRB nature. Using data from Einstein, 



ROSAT and ASCA, 



Larson & Schulman 



days. This was confirmed by 



Dubusetal 



1997) have suggested that the actual period is 3.4531 



1999), who reported a period of 3.4535 +0.0005 days. 
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Pietsch et al 



(120041) have analyzed X-ray data from XMM-Newton and Chandra and derived 
an improved period of 3.45376 + 0.00021 days. They also used optical B and V data to suggest 
the optical companion is of spectral type between B0 I and 07 I, with masses of 25-35 Mq, 
making the compact binary component a black hole. If this is true, X-7 is the first eclipsing black 
hole high mass X-ray binary (HMXB). 



In the course of the Chandra ACIS-I survey of M33 (ChASeM33, 



ray eclipse ingress and egress of X-7 were observed for the first time ( 



Sasaki 


2005 


), the X 


( 


Sasaki et al. 


2005) 



Pietsch et al 



(120061) combined Einstein, ROSAT and XMM-Newton data, and derived an im- 
proved period, of 3.453014 ± 0.000 020 days and a possible period decay rate of P/P 
-4 x 10~ 6 yr -1 . 



Pietsch et al 



association HS 13 



2006) have also used HST WFPC2 images to resolve the OB 



Humphreys & Sandage 



1980» in which X-7 resides, and have identified the 



optical counterpart to be an 06 III star whose mass is at least 20 Mq. Based on the optical coun- 
terpart identification, possible o rbital period decay r a te, lack of pu l sation s, X-ray spectrum and 



Pietsch et al 



(2004). 



Pietsch et al 



(2006) suggested that X-7 is 



fits to the optical light curves of [j 
a black hole HMXB. 

In the previous studies the periodicity of X-7 was derived from the X-ray data, and no in- 
dependent analysis o f optical data was pe rformed. Optical light curves of X-7 were presented 



in B and V filters by 



Pietsch et al 



(2004), who folded the optical measurements on the X-ray 



pe riod. This paper prese nts an analysis of optical light curves of the X-7 counterpart obtained 



by 



Hartman et al 



(120061) in a search for variable stars in M33, performed at the C anada-France 



Hawaii-Telescope (CFHT) with the Sloan r' and i' filters 



are combined here with the B and V observations of Pietsch et al 



Hartman et al 



(2006) measurements 



(2004) to search independently 



for a periodic modulation. The combined data, as presented in §|2]do show clear pe riodicity with 
the X- r ay period, as presen ted i n § [3] Based on broa d-band UBVRr'i' data from 



(2006), 



Massev et al 



d2006l) and 



Pietsch et al 



Hartman et al 



( 120061) we apply a photometric modelling to de- 



rive the stellar temperature and radius of the optical component in § 14.11 and periodic ellipsoidal 
model in § 14.21 We briefly discuss our results in §|5] 



2. The Optical data 

The present analysis is based on the optical measurements of X-7 in four bands, obtained from 



Pietsch et al 



two se ts of observations carried out ~4 years apart. The B and V light curve s of 
( 2004 , Table 4) were obtained from a special photometric analysis of DIRECT^ ( Mocheiska et al 
200 lb images, taken at the Kitt Peak National Observatory 2.1-m telescope in 1999. The V light 
curve includes 70 measurements and the B light curve only 30 points. Sloan r' an d i' light curves 



ve\PI Jllai 



20061 object 243718). Their data were obtained at the CFHT 3.6-m telescope atop Mauna Kea on 



two consecutive observing seasons in 2003 and 2004. An i' finding chart of X-7 is given in Fig.Q] 



1 http://cfa-www.harvard.edu/~kstanek/DIRECT/ 

2 http://www.astro.livjm.ac.uk/~dfb/M33/ 
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34.4 1:33:34.0 33.6 
Right Ascension (J2000) 



Fig. 1. X-7 finding chart, in i', from the CFHT data. Position of the variable object is indicated 
by an arrow. Image FOV is 20 x 20 arcsec. North is up and East is to the left. 



A comparison with the finding chart of 



Pietsch et al 



2006 , their figure 5) shows this object is at 



Pietsch et al 



(2006). The identification of the 



the position of X-7 optical counterpart identified by ] 
same period in the optical as in the X-rays (presented hereafter) confirms the identification. The 
i' light curve contains 32 measurements, where we ignored two measurements with uncertainty 
larger than 0.04 mag. The r' light curve consists of 31 measurements, where we ignored three 
obvious outliers. 



3. Period analysis 



To search for period icity we have applied a multi-band double-harmonic fitting 



(iShporer & Mazehl 120061) to the four available light curves together, which included 158 
individual measurements. For each trial period we fitted the entire data set with a double- 
harmonic function, with independent zero points for each of the four light curves. The 
periodogram value for each trial period was taken as the power (sum of squared harmonic 
coefficients) divided by the x 2 goodness-of-fit parameter of the fitted light curve. 

Fig. represents the resulting periodogram. The strongest peak corresponds to a period of P opt 
= 3.4530 + 0.0014 days . The second strongest peak is exactly at the first harmonic of P opt . The 
third strongest peak is at a frequency corresponding to the lunar cycle. T he optical per i od ag rees 



Pietsch et al 



(2006). 



well with the known X-rays period, P x = 3.453014 + 0.000020 days, of|] 

Fig. [3] presents the four light curves folded o n the detected period , with phase zero taken at 
HJD = 2453 639.119 — the mid- eclipse time of ; 
and V data of 1 



Pietsch et al 



Pietsch et al 



(2006). The figure shows the B 



(120041) for compl eteness, including seven V measurements that were 



Pietsch et al 



2004 since that figure brings only V measurements that 



not plotted in figure 4 of 
were obtained close to B measurements. Fig. [3] clearly shows ellipsoidal modulation in all four 
bands. For each band independently, we fitted the data with two harmonics and derived the total 
amplitude of the modulation and the phase of minimum. The results are listed in TableQ] For the 
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Fig. 2. Multi-band double-harmonic periodogram. 

Table 1. Harmonic fits amplitude and minimum phase near 0.0. 



Filter 


Amp. 
[mmag] 




Min. Phase 




i' 


47 + 3 (63 


±4) 


0.962 ± 0.006 (0.949 


+ 0.006) 


r' 


39 ± 2 (48 


±2) 


0.960 ± 0.003 (0.953 


± 0.003) 


V 


34+1 




0.970 ± 0.001 




B 


30+1 




0.979 ± 0.001 





i' & r' bands we ignored the four faintest points in each of these bands. The amplitudes derived 
when we included these points in the analysis are given in parenthesis. 

It appears as if the amplitude of the modulation decreases with decreasing wavelength, from i' 
to B, even when ignoring the four faintest points in ;' & r' . The origin of this varying modulation 
is not clear. It could result from some real wavelength dependence of the system's modulation. 
On the other hand, the different modulation could come either from some secular changes during 
the -4 years between the KPNO B & V and CFHT /' & r' observations, or from some wavelength- 
dependent blending of light from a near-by star. 

4. Modelling the system 

In our modelling of the X-7 binary system, we preferred to find first a stellar model that fits the 
broad-band photometry of the system, assuming the optical brightness is coming from the optical 
star alone. From the stellar radius, the observed X-ray eclipse width and the periodic ellipsoidal 
modulation we can derive the two masses as a function of the orbital inclination. 
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Fig. 3. Folded light curves, in magnitude, in (top to bottom) ;', r', V and B. All light curves are 
folded on the period identified here, of 3.4530 days, and phase zero was tak en to be at HJD = 



2453639.1 19, which is time of mid-eclipse in the X-rays (Piets ch et al.1 12006) 



4.1. Optical star modelling 



We h ave used the following broad-band magnitudes: U = 18 .1, B - 18.8, V = 18.9 (IPietsch et al 



2006, based on HST data ), R = 19.0 dMassey et al 



2006, LGG 



mag dHartman et al 



A r' 



19.1 and i' = 19.42 



2006), assuming an uncertainty of 0.1 mag on all values. We converted the 



observed measuremen ts to absolute magnitudes by using an M33 distanc e modulus of 24.62 mag 



( Freedman et 



Schlegel et al 



(1998). 



2001), A v = 0.53 + 0.06 mag from lPietsch et al.1 00061) and A A /A V values from 



We compared these m agnitudes to the ones determined fr om the integrated m agnitudes from 



Girardi et al 



(2002. 



2004), where we assume M bo i Q = 4.77 dGirardi et al 



2002) and fit for the 



radius for each value of log(g) and T e jj. We did this for both [Fe/H] = and [Fe/H] = -0.5 tables, 



3 The Local Group Galaxy Survey, http://www.lowell.edu/users/massey/lgsurvey/index.html 
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Fig. 4. Photometric broad-band modelling of the X-7 stellar component. Error bars present the 
assumed uncertainty of 0.1 mag in all bands. Solid line presents a model with T e ff = 27000 K, 
logg = 5.0 and R = 25 Rq. The dashed and dotted lines the models at 95% confidence boundary 
of our modelling. The dotted line is a model with T e ff = 18000 K, R — 35 Rq, the dashed line is 
a model with T eff = 47000 K, R = 15 Rq. 



since assuming X-7 is at a galactocentric distance of 2 kpc results in [Fe/H] = -0. 22, based on 



Garnett et al 



d 19971) and the [O/H] to [Fe/H] conversion of 



Maciel et al 



the [O /H] gradient of 
J20031 

The data allow a large range of stellar models. One such model, for [Fe/H] = 0, with x 2 = 
1.62, has T eff = 26000 K, logg = 5.0 and/? = 25 Rq. For [Fe/H] = -0.5, with X 2 = 1-55, we get 
T e ff = 27000 K, log g — 5.0 and R — 25 Rq. The 95% uncertainty range, where x 1 increases by 
up to 8.0 (for three parameters), is: 18000 < T eff < 47000 K, 15 < R < 35 Rq. The value of log g 
is unconstrained. Fig.|4]brings three representative models; the solid line presents the model with 
T ef f = 27000 K, R = 25 Rq, the dotted line represent the model with T eff = 18000 K, R = 35 Rq 
and the dashed line shows the model with T e ff = 47000 K, R — 15 Rq. The figure shows that the 
absolute magnitudes of the X-7 optical counterpart can not constrain the stellar radius, and any 
value in the range of 15-35 Rq is possible. However, as we show in the next section, a stellar 
radius of R > 20 Rq is very unlikely, since the star will be too massive. Therefore, our most 



likely models are in the range of 15 < R < 20 Rq and 33000 < T eff < 47000 K. 



4.2. Light-curve modelling 



We used two different codes, PHOEBE and ELC, to model the obtained periodic light curves and 
derive an estimate of the masses of both components. For both codes, the main effect is due to 
the tidally induced ellipsoidal shape of the optical component. 

The first program is the "PHysics Of Eclipsing Bi nariEs" program of |Prsa & Zwitterl (120051) 



a fron t-end code for the Wilson-Devi nney program dWilson & Devinnev 



19901) . Following 



Pietsch et al 



1971 



Wilson 



1979 



(12006b . who used PHOEBE to model the X-7 B and V light 
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curves, we used the eclipse half angle value of 26.5° + 1.1° to find the ratio of the radius of 
the optical component to the semi-major axis as a function of inclination. We varied the incli- 
nation between 90° and about 70° and fit for the mass ratio and a luminosity scale to match 
each light curve, as suming a temperature of 27000 K and linear limb-darkening coefficients from 



Claret (2000 



20041) . assuming a metallicity of [Fe/H] = -0.2. 



We performed the fits ignoring the four faintest points in each of the r' and i' light curves 
which occur near phases and 0.5 — no model was able to match these points. Including those 
points does not substantially change the mass-ratio, while it does slightly affect the luminosity 
scale. 

There is not enough data to constrain the stellar radius or the inclination of the system; a good 
fit to the light curves could be found for large ranges of these two parameters. One such model is 
plotted in Fig. [5] The PHOEBE light curve is drawn with a dotted line and the other model with 
a dashed line. For both, values of R — 20 Rq, T e ff = 33000 K and i = 80° were taken, although 
changing these values results in essentially no difference for the model light curve. 

The range of allowed inclinations and stellar radii did not enable us to constrain the masses 
of the two components and they could only be derived as a function of the assumed orbital 
inclination and optical radius. We found that for the PHOEBE code the mass ratio (M x /M opt ) 
varies from 0.29 at 90° to 0. 15 at 70°. Below 70° the optical component would exceed its Roche 
lobe. Assuming R = 20 Rq, the mass of the optical component would vary between 79 Mq and 
49 Mq, while that of the compact object would vary between 23 Mq and 7.5 Mq. Because the 
fits are fairly insensitive to the temperature of the optical component, the resulting masses scale 
as the radius of the star to the third power. For R = 15 Rq, the mass of the optical component 
varies between 33 Mq and 21 Mq, while the compact object mass varies between 10 Mq and 
3 Mq. We note, however, that assuming a stellar radius larger than 20 Rq results in an unrealistic 
large stellar mass. We therefore limit the discussi on to stellar radii sma l ler tha n 20 Rq. 

The second program used was the ELC code o flOrosz & Hauschildtl (120001) . This program has 
been specifically tailored to model optical light curves of X-ray binary systems, including effects 
such as the presence of an accretion disk, X-ray heating of the optical component and the ability 
to explicitly take one of the components to be invisible. To model the X-ray heating we assumed 
an average X-ray lum inosity of 1Q- 37 5 erg s _1 , and to model the disk we assumed an inner disk 
temperature of 10 7 K ( Pietsch et alj|2006 ). We assumed a disk opening angle of 2°, and a power- 
law exponent for the temperature profile of the disk of -0.75. We used the black-body mode of 
the model rather than using the detailed model atmosphere since the model atmosphere would 
not necessarily be a good match to the Sloan r' or ;' filters. We proceeded as above, stepping 
through a range of inclination angles for two values of the optical radius and fitting for the mass 
ratio. This time we also fit for the inner and outer disk radii. One possible model, for R — 20 Rq, 
T eff = 33000 K and i = 80°, is plotted in Fig.S 

The resulting mass ratio varies from 0.23 at 90° inclination to 0.07 at 65° inclination. For a 
radius of 20 Rq this corresponds to a mass range of 83 to 42 Mq for the optical component and 
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Pietsch et al 



Fig. 5. BVr'i' light curves phased at P = 3.453014 day and T = 2453639. 1 19 from 
I2OO6J), with model light curves for R = 20 Rq, T eff = 33000 K and i = 80°, overplotted. The 
dotted line is the PHOEBE light curve, with no accretion disk or X-ray heating. The dashed line 
is the ELC model which includes an accretion disk and X-ray heating. 



a range of 19 to 3 Mq for the compact object. Below 65° the ELC model star would exceed its 
Roche lobe. In units of the inner Lagrange radius of the compact object, the inner disk radius 
varies from 1.6 x 10~ 5 to 5.4 x 10~ 5 while the outer disk radius varies from 0.61 to 0.81. 

The masses of both components as a function of inclination, derived by both codes for two 
values of the optical radius, are plotted in Fig. [6] 



5. Discussion 

We have conducted a broad-band optical analysis of the X-7 optical counterpart, including period 
analysis, photometric modelling and light curve modelling. We applied a multi-band double- 
harmonic period analysis of X-7 optical light curves in four bands and identified a period of 
P op t = 3.4530 + 0.0014 days . This period, de rived from an indep endent optical analysis, is in 



good agreement with the known X-rays period dPietsch et al . 2006). 



Using photometric absolute magnitudes in six bands, we find a range of models for the optical 
counterpart, with 15 < R < 35 Rq and 18000 < T ef f < 47000 K. However, stellar radius larger 
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Fig. 6. Plots showing the mass vs. inclination for two values of the stellar radius. The left plot is 
for the optical component, the right plot is for the compact object. The solid and the dashed lines 
show the results for R — 20 Rq, using the ELC and the PHOEBE codes, respectively. The dotted 
and the dot-dashed lines show the results for R = 15 Rq, using the ELC and the PHOEBE codes, 
respectively. 



than 20 Rq results in an unrealistic large stellar mass. We therefore suggest that the stellar radius 
is smaller than 20 Rq. The possible temperature range is reduced accordingly to 33000 < T e ff < 



Pietsch et al 



(2006) is well within our 



47000. The classification of the optical star as 06 III by 
uncertainties. 

We have modelled the optical light curve using two programs. While the ELC model, incor- 
porating both a disk and X-ray heating matches better to the data than the PHOEBE model (see 
Fig. |5j, neither model yielded a good fit to the light curves. In particular, the minima in r' and i' 
appear to be underestimated by all models that we have tried. At this point it is unclear if this is 
due to random or systematic errors in the observations or to inadequacies in the physical models. 

The simplistic models can nevertheless yield a rough estimate of the masses of the two com- 
ponents as a function of the optical radius and the orbital inclination. For the smallest likely 
optical radius (15 Rq) we get 3 Mq and 1.3 Mq from the PHOEBE and ELC codes respectively. 
Therefore, the present analysis can not exclude completely a neutron star as the compact object, 
although this option is very unlikely. 

On the other hand, the similarity between M33 X-7 and the three known black-hole HM XB 
— Cyg X-l, LMC X-l and LMC X-3 (e.g., 



Cowlev 



1992 



McClintock & Remillard 



2003). is 



striking. All four systems have orbital periods of a few days, their optical counterpart is an early- 
type star and the compact object has a mass of 6-10 Mq. Radial velocity measurements of M33 
X-7 will enable us to better constrain the masses of the two components. Large telescopes and 
presently efficient spectrographs render such measurements feasible. Together with the X-ray 
eclipse and the optical light curves, we should be able to understand this system better than any 
other black-hole HMXB. 
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Finally, if this system would be proven to be a black-hole HMXB, it is interesting to note that 
three out of the four such known systems were found in external galaxies. This might indicate 
that we are not very efficient in detecting such systems in our own galaxy. Maybe the dust in the 
Galactic plane hides from our telescopes many more Galactic black holes. 
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